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Abstract: We consider a model of inflation which has recently been proposed in the 
literature and where inflation is induced by corrections to the energy density coming from 
the non-commutativity of spacetime. We show that the very rapid inflationary expansion 
typical of this model is responsible for a burst of particle production which ends inflation 
and leads to a radiation-dominated phase. We analytically estimate the energy density of 
these particles and we confront the results with more precise numerical calculations. We 
estimate the number of inflationary e-folds before the back-reaction of the radiation energy 
density overcomes the non-commutative effects and terminates inflation naturally. 
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1 Introduction 

The inflationary paradigm is based on the assumption that the early Universe has experi- 
enced a phase of quasi-exponential expansion, which is usually driven by the dynamics of 
a minimally coupled scalar field [1-3]. The potential associated to the scalar field and/or 
the initial condition must be carefully modeled so that the inflationary expansion lasts suf- 
ficiently long to solve the horizon and flatness problems and, at the end of inflation, the 
scalar field can decay into a sufficient amount of matter. Even though simple single field 
inflation is in good agreement with observations [4], it remains a phenomenological model 
as the origin of the scalar field is unknown. Another serious problem is the large degree of 
extrapolation of physical laws, that are extended to the extreme, near-Planckian regime at 
which many models of inflation take place. 

It is a well-known fact that General Relativity (GR), as it stands, no longer holds 
when the energy density reaches Planckian values. Several attempts to unify GR and 
quantum mechanics, in order to find a consistent quantum theory of gravity, are under 
investigation. Among these, there is one line of research that suggests the extension of 
the Heisenberg uncertainty principle to pairs of coordinates. In this theory, the quantum 
algebra relating the position and the momentum of a single particle is implemented by a 
new non-commutative relation of the form = iOij, where 9ij is of the order of the 

Planck length squared [5]. This commutator clearly restricts non-commutative effects to 
very special situations, where the energy density is of the order of the Planck value. Such 
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extreme conditions can occur near singularities and at the onset of an inflationary phase. 
Non-commutative effects in these situations have been already investigated, see e. g. [6]-[14]. 

In this paper we consider a model of non-commutative inflation proposed in Ref. [15], 
which is based on the so-called coherent state formalism developed in Refs. [16]. The model 
predicts a rapid phase on inflation driven uniquely by non commutative effects, that act 
on the quantum spacetime fluctuations, modeled as a scalar field, in such a way that they 
behave as a time-dependent dark energy term with a Gaussian density profile. 

The advantage of this model is that it does not require any ad hoc potential and the 
inflationary phase ends naturally when the energy density drops below a certain value. In 
the simplest setup, any form of pre-existing matter is absent and the Universe is initially 
empty at t — )• — oo. In this paper we study particle creation in this time-dependent back- 
ground at high curvature. The principle of this effect was studied long ago, and it is known 
to be quite effective in a quasi-de Sitter space [17]- [20]. In the standard inflationary sce- 
nario, the production of hot and relativistic particles is subdominant during inflation and 
it occurs only after the slow roll phase, when the inflaton starts oscillating, at the onset of 
a re-heating or pre-heating phase. 

In our model there is no inflaton therefore quantum particle production is the only way 
to fill the Universe with radiation. In fact, as we will show, the non-commutative expansion 
is so rapid that particle creation is actually "explosive" and it eventually ends inflation by its 
back-reaction on the spacetime. In this paper we study this process in detail and calculate 
the energy density of the generated particles. 

In Sec. 2 we briefly review the non-commutative inflationary model at hand and in 
Sec. 3 we give an analytic estimate of the particle production. In Sec. 4 we present the 
numerical results, and we conclude in Sec. 5 with some remarks. 

2 Non-commutative inflation 

In Ref. [15] one of us has shown that, in a homogeneous and isotropic Universe, non- 
commutative effects lead to a dynamical vacuum energy density that has the form 

(p)Nc-/Ooexp(-AtV^), (2.1) 

where 6 has the dimension of a length squared. In static spacetimes, a similar formula holds 
for the energy density for the mass at the center of a black hole, see [21]. As shown in [15], 
these results follow from a redefinition of relativistic quantum fields and the value of 9 turns 
out to be of the order of the Planck length squared £p. Remarkably, the last statement is 
consistent with traditional vacuum energy calculations. In fact, a simple argument shows 
that the expectation value for the vacuum energy can be calculated with the integral 

= i ^^(2^2^'^^^' ^^-^^ 

which yields, in the limit A 3> m, (p) = A^/(167r^). The ultraviolet cutoff A identifies the 
transplanckian frontier and it can be set to A = (SvrG)"^/^, where G is the Newton constant 
(with = c = 1), see e.g. [22, 23]. 
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In the coherent state non-commutative approach, the vacuum energy for the fields 
should be computed instead as [9] 



(p)nc=/ dk -Vk^ + m\ (2.3) 

Jo (271-)^ 2 

where now 9 is the non-commutative parameter that corresponds to a length squared. In 
the limit where m is negligible, this integral gives (/5)nc = 1/(80^). By equating the two 
results, one finds that 

Ve = (128)^/V£p ~ 10.6 ^p, (2.4) 

where £p = \fG is the Planck length. The two calculations are then consistent provided 
the parameter \fQ is of the order of the Planck length. To account for this result, we will 
parameterize Q by 

\/^ = a^p . (2.5) 

Naturally, we expect a to be a constant of order unity. In the derivation given in Ref. [15] 
it was found that the leading term of the deWitt-Schwinger expansion of the energy mo- 
mentum tensor for a single massless scalar field in the non-commutative field theory, based 
on the coherent state formalism, is 

(^/Lti/) — 2271-2^2 5/^*" ' 

which corresponds to an energy density somewhat smaller than the above estimation. How- 
ever, there is the possibility that the vacuum is determined by several fields increasing the 
value of its energy density. 

The formula (2.1) yields the Friedmann equation 

= ^Pnc(t) = i^^e-(*-*-)^/^\ (2.7) 

where a is the scale factor of the metric As^ = —dt^ + a?'{t)5ijdx'^dx^ and H = a' /a. From 
now on, the prime indicates the derivative with respect to the physical time t while a dot 
with respect to the conformal time ry, defined by dt = adrj. Without loss of generality, 
we can also set tmax = 0. The parameter Hm absorbs other eventual factors, such as 
the number of fundamental fields that are excited by non-commutative effects, and it is 
time-independent. If we equate the niaxinium energy density, Pmax — 

p{0), with the (r^) 

component of the expression (2.6) we find 

In the following we shall, however, relax this condition on Hm and let /3 be arbitrary. As 
mentioned above, a larger value for /3 is easily justified by allowing for N massless scalar 
fields contribution to Praaxi 

such that N/]"^. 
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Like in the black hole case, Eq. (2.7) has to be interpreted as the impossibility of 
localizing an arbitrary amount of energy at the time t = as the result, for example, of a 
collapse. However, to simplify our calculations, we choose initial conditions supposing that 
the Universe starts in the far past empty and with a small value of a and of the energy 
density. Thus, we only consider the non-commutative energy density and neglect any other 
contribution from radiation or cosmological constant. This is one of the advantages of this 
model: we can let the system evolve from very simple vacuum initial conditions. The only 
arbitrary choice is the sign of H, which we take to be positive. 

With these hypothesis, we easily find the solution to Eq. (2.7), which reads 



a{t) = am exp 



Krr, V27:e Erf 



t 



2V2e 



am exp 



27r/3a Erf 



t 



2V2e 



(2.9) 



Erf(x') = / e-^'dx , 
Jo 



where am is an integration constant. It can be checked that a" is initially positive and then 
changes sign at the time when the comoving Hubble length (aH)~^ reaches its minimum, 
while H reaches its maximum, Hm, at t = 0, see fig. 1. In other words, the global evo- 
lution of the scale factor is very similar to the case of the pre-Big Bang scenario of string 
cosmology [24]. H'(t) = —^H{t) is positive for t < and negative for t > 0. 
The number of e-foldings of inflation can be estimated as 



iVe ^ In f ^ = JSttHIO = Vs^a^ . (2.10) 

\a[t — )• — oo) / 



If we insert /3 from (2.8), requiring 60 e-folds as in standard inflation, yields a ~ 0.014. 
Leaving both, a and (3 free, the requirement of at least 60 e-folds of inflation implies the 
relation 

60 

a/3>^^~12. (2.11) 
x/8tt 



Note that this is the 'naively estimated' number of e-folds. In Sec. 4 we shall show that the 
true number of e-folds is reduced by up to a factor of 2 due to backreaction from particle 
creation. 

By imposing the second Friedmann equation, p' + 'iH{p -|-p) = 0, we find the effective 
pressure p. Defining the equation of state parameter a; by p{t) = u;{t)p{t), we obtain 

This leads to a dynamical crossing of the value oj = —1 value at t = where H{t) reaches 
its maximum. At negative times, this models exhibits an effective "phantom" behavior, 
w < -1. 

3 Particle production 

It is a long-known fact that a sufficiently rapid growth of the scale factor leads to particle 
production [17]- [20]. Technically, a meaningful calculation requires two asymptotic regimes 
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Figure 1. Qualitative behavior of a (solid line), a" (dotted line), H (dashed line), and (aiJ) 
(dot-dashed line) as functions of time. 



where the the spacetime is flat, or slowly varying, so that the particle number is a well- 
defined quantity. In our simple model the scale factor is asymptotically constant both in 
the past and in the future, an ideal situation to calculate the energy density of the produced 
particles. 

The mode equation cannot be solved analytically because of the complicated form of 
the scale factor (2.9), so we need to resort to numerical approximations. However, some 
analytic estimates can be given. 

Let us consider the action for a massless scalar field 

S = -\j d''xV^g^''d^^d,<p. (3.1) 

On a Friedmann-Lemaitre (FL) background with metric ds'^ = a{'rff'{—dif' + dx^) in con- 
formal time, it becomes 



S = d xdr] 



2 

Setting xiv) — this action takes the canonical form 



(3.2) 



S = lj d'lLdn [x' - (Vx)' + Ix'^j , (3.3) 
and the associated equation of motion, in Fourier space, reads 

Xk + ^kiv)Xk = 0, where a;|(r?) = k'^ - d/a = k'^ + m'^{ri) . (3.4) 
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At early times Tn'^{r]) < indicating an instability and at late times it becomes positive. 
Following the usual quantization procedure, the field x is promoted to an operator Xkiv) — 
(^kvlirj) + a^_i^Vk{i]) and the modes v satisfy the classical equation of motion 



Vk +Lol{r])vk = 0. 



(3.5) 



In our case, we find it more convenient to work in physical time, since the form of the scale 
factor and its derivatives are given explicitely in t. The above equation then takes the form 



v"+[-\v' + 
a 



i\ 2 



u = 0. 



In terms of the dimensionless time variable T = t/\/89, momentum K = k 
dimensionless mode function V = v/{86)^^'^ the mode equation finally reads 



(3.6) 



, and the 



d^V 1 da dV 
dJ^ ^ adTdT ^ 



K 



1 f d^a 



dT^ 



1 / da\- 
^ ydTj 



V = 0, 



(3.7) 



and the initial vacuum state, corresponding to a = constant, is ^ = exp{—iKT)/^/2K. In 
the variable T, the scale factor and the background energy density are given by 



a(r) = a^exp 



Erf(r) 



4o2 



Pnc(r) = Pmaxexp(-2r^), Pi, 



3 a*/3 



3 N^a^ 
647r2 02 



(3.8) 
(3.9) 



We also report the expression of the dimensionless dynamical mass M^{rj) defined in Eq. 
(3.4), namely 



TT 



exp 



Ne Erf(r) 



2 7^ exp(r J 



(3.10) 



3.1 Definitions of the energy density 

The usual definition of the energy-momentum tensor in GR is 

2 S{y^\CM) 



5g 



(3.11) 



where C m is the matter Lagrangian density. This is the definition of the energy-momentum 
density that has to be used in the Einstein equations. In general, it does not agree with the 
canonical energy-momentum tensor obtained via the Noether theorem, but differs from it 
only by a total derivative. 

The Hamiltonian related to the action (3.2) is 



(3.12) 
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In the previous sections, we found useful to change the variable (/>—)• x = CL(p, which is a 
canonical transformation, in order to obtain the standard equation of motion (3.4). The 
action for x given in eq. (3.3) leads to the Hamiltonian 

= ^ / rf'x (x^ + (Vx)^ - V) . (3.13) 

differs from since the canonical transformation — )• x is time-dependent. The 
difference between the two Hamiltonians is given by the time derivative of the generat- 
ing function of the canonical transformation [25]. In any case, must be used to find 
the energy density on the right hand side of the Friedmann equation because this is the 
same energy associated to the energy-momentum tensor obtained by varying the matter 
Lagrangian with respect to the metric g^,^. Nevertheless, out of interest, we shall also 
compute the energy density associated to H^. in our numerical analysis. 

3.1.1 Adiabatic subtraction 

In the previous paragraph we have defined the classical expression of the energy-momentum 
tensor T^,^. After quantization, one obtains the corresponding quantum operator T^,^. The 
semi-classical approach consists in replacing the right-hand side of Einstein equation by 
the expectation value {s \ T^u \ s) where | s) is the quantum state of the scalar field. 
Without the exponential ultraviolet cutoff coming from non-commutativity, see Eq. (2.3), 
this expectation value is infinite and needs to be regularized. For example we have 

(0 I HM I 0) = ^5(3)(o) j [\vk{v)? + u^lmvkiri)?] . (3.14) 

This is the definition of energy density discussed in Ref. [29]. The spatial energy density is 
obtained by dividing this expression by the 'volume' (27r)^a^(5(^^(0). Performing the angular 
integral, the unrenormalized energy density becomes 

/^r = 4^ / dkk^[\v,{r^)\^+^l{r^)\vM\'']- (3-15) 

From this quantity we have to subtract the (flat-space) vacuum contribution which is the 
same integral as above but with Vk equal to the plane- wave expression, Vk = ex.p{—ikri) j \plk. 
This "zero-point subtraction" is not sufficient when the background is curved as the vacuum- 
subtracted expression, p^<^ is still divergent and one needs to resort to more sophisticated 
techniques, see e. g. [26]. There are several proposals in the literature on how to apply 
this regularization and on its domain of validity. In particular, adiabatic subtraction has 
been recently discussed in relation to the regularization of the primordial power spectrum 
[27, 28]. Usually, one subtracts from the integrand of the divergent energy density its ex- 
pansion in WKB, slowly varying, terms. The lowest order WKB term coincides with the 
"instantaneous vacuum solution". The subtraction needs to be performed up to the 4th 
WKB order to have a finite result in the ultraviolet (UV) regime. 

However, when the time dependence of the background is very strong, the different 
subtraction schemes give very different, sometimes even negative results for the energy 
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density and they cannot be trusted. In fact, when the time variation of the background is 
very rapid, the notions of particle and energy are simply not well defined. In our model, 
all the definitions proposed in the literature give the same result at late times, when the 
background is very slowly expanding, and we can reliably compute the energy density at 
the end of inflation. 

For the field (/>, the energy density upon adiabatic subtraction in conformal time is (see 
appendix A) 



p<i>{v) 



1 

47rV 



j dk 


a2 




2 












1 


aa^ - 






8P 


a3 



+ 



A?5 



(3.16) 



where the last three terms in the square bracket correspond to the counter-terms of adiabatic 
order (the term —k), 2 (the term oc 1/k) and 4 (the term oc l/A;^). Both the 0th and 2nd 
order terms dominate in the UV, where we expect adiabaticity to be a valid approximation, 
so that there is a reasonable justification to subtract these terms. However, the 4th order 
term is logarithmically divergent, hence it dominates also in the infrared (IR) regime, where 
adiabaticity is not at all verified. In Ref. [27], the subtraction is justified with the fact that 
the final result converges. However, we still believe that it is not a good expression for 
the energy density of the particles generated during inflation [28]. In particular, the 4th 
order subtraction introduces a new infrared divergence which is unphysical and which is not 
present in the un-subtracted result. Furthermore, as we shall see in our numerical results, 
this energy density can become negative, a most unphysical behavior for the energy density 
of a collection of free relativistic particles. 
In physical time, Eq. (3.16) becomes 



1 



k - ia' 



1 dk k^ 


a' 











12 



1 



2k 8A;3 



/ 2 //2 I o /4 r, 2 III I n 12 ll\ 

(a a + 3a — 2a a a — 2a a ) 



(3.17) 



The notation /)0^„ indicates that we have performed the adiabatic subtraction up to 
the order n. By repeating the same calculations for Hy., we find, in conformal time 



Pxiv) 



1 



47r2a4 



dk k^ 



^ 1 a 1 

2ka ^W\a 



\vk\'+ik'-- 
a 



\Vk\^ + 



(3.18) 



+ 



k^ 



Again, the three last terms in the square brackets correspond to the adiabatic counter-terms 
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of order 0, 2 and 4. In physical time the above expression becomes 

1 



47r2a4 



dkk' 



\av 



.'f + {k^-{a'' + aa"))\vkf + 



7 I ^ / i2 , //\ I ^ / /2 I //\2 



(3.19) 



As explained above, these definitions make sense only if the time dependence is weak, 
i.e. when the spectrum is dominated by frequencies such that k ^ a/a. In this situation, 
the adiabatic counter-terms are small and all the expressions for the energy density do not 
differ by much. In opposition, in the infrared regime the adiabatic approximation breaks 
down and the result cannot be trusted. 

In the context of the non-commutative regularization, the UV divergences are cured 
by construction. Indeed, the exponential damping term in the integral measure ensures 
the convergence of the energy density for /c — )• oo. Thus, to be fully consistent with the 
underlying non-commutative field theory, the above expressions for the energy density have 
to be modified in the UV where the exponential damping from non-commutativity sets in, 
see [15]. More precisely, eqs. (3.17) and (3.19) have to be replaced by 



1 



47r2a4 



- k 



12 



dk k exp 
1 
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[a a + 6a 



2k 



Sk^ 



/4 



n 2 III I 

la a a 



n l2 II \ 

2a a 



(3.20) 



and 



47r2a4 



j dk k^ exp 









|af^|^ + {k^ ~ [o!"^ + aa")) \vk\ 



; I ^ t i2 , II \ , ^ / /2 I II \2 

k-\-—ia + aa ) + —ri^ia + aa 



2k 



8k^ 



(3.21) 



3.1.2 Analytical estimation 

Before presenting our numerical results, we perform an analytic approximation for the 
energy p^. The un-subtracted, ultraviolet finite energy density is 



1 



2^4 



y dk k^ exp 






a^ 




2 9, ,9" 

+ k^\vk\ 













(3.22) 



Let us consider a mode with wave-number k that originates as a sub-Hubble scale in the 
far past, crosses the Hubble horizon at rjo^tik) and finally re-enters at r][^{k), see fig. 2. 

As long as the mode is well inside the horizon, we have uj'^ = k'^—d/a ~ fc2, its amplitude 
can be approximated by a constant and the energy density can be set to zero. After 
crossing the horizon, the mode starts to grow proportionally to the scale factor according 
to (vk/a)' = 0. To show this, note that in the super-horizon regime we have a;2 ~ —a/a 
and the equation of motion (3.5) reduces to 



Vk 



a 
a 



(3.23) 
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/ y aa " + a' ^ 1 












Tatit 


T 

in 



Figure 2. The dynamical mass \m^{T)\ (blue line) and a co-moving wave number (red line) are 
indicated as functions of the dimensionless physical time T. When k enters the shaded region the 
corresponding wavelength exits the Hubble scale, at Tout < 0, and it re-enters it when k crosses 



the Hubble scale again at Tin > 0. At T 
particle creation is halted for all modes. 



1, 



becomes positive and the instability leading to 



which has the general solution 

Vkiv) = ci{k)a{r]) + C2{k)a{r]) 



dT]' 



VO 



a?{rj) 



(3.24) 



Since the second term decays rapidly, we soon obtain Vk oc o. In the computation of the 
energy density this leads to the approximation 



k a{r]y 



(3.25) 



that holds when the mode is super-horizon, namely for rjontik) < rj < r/in(A;). When 
V > Vinik), the mode is again sub-horizon and thus of constant amplitude. With this 
approximation, all modes with wave number above /Cmax = Y^max(a/a) are never excited 
and the energy density is always finite. We then obtain the following approximation for p^: 



Pest 



47r2a4 / 



dk k exp 



-26 



k 



(3.26) 



with 



fik,v) 



1 



V < Vout{k) 
a(ry)/o(r/out(fc)) r]out{k) < 7] < Vm{k) 
aiVinik)) / a{riont{k)) V > r/in(fc) 



The term —k/2 in the integrand comes from the requirement that the energy density is zero 
in the far past and that it is continuous. 

In the limit a — t- constant, the definitions and p^ are the same and coincide with 
energy density of the field in Minkowski space. This is can be traced back to the fact that 
for a =constant, the canonical transformation is time-independent so the two Hamiltonian 
densities are the same. In the opposite regime, when T-L^,T-l > k"^, the adiabatic subtraction 
is no longer valid and and the counter-terms diverge, causing to become most likely 
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negative. For p-^ the situation is even worse: the Hamiltonian H^^ has no minimum when 
< which corresponds to the non-adiabatic regime [29]. Negative energy densities are 
inevitable, with the formulae above, because a consistent definition of particle and energy 
density is not known in strongly time-dependent backgrounds. However, as we shall see 
below, these expressions all coincide after the inflationary phase, when H ^ 0. Hence they 
can give a reasonable estimate of the final energy density of scalar particles. 

As soon as becomes of the order of the non-commutative energy density (3.9), back- 
reaction becomes important and we can no longer trust our expansion law since the contri- 
bution from Pfjj in the Friedmann equation cannot be neglected. The late non-commutative 
phase is altered by the energy of the particles produced in the early phase in such a way 
that the overall, final energy density can be different from the estimated one. 

An advantage of the approximation pest described above is that it is always positive and 
finite. However it may not be very precise, mainly because the modes do not truly switch 
from the sub-horizon to the super-horizon regime instantaneously. In the non-commutative 
model at hand, the approximation is expected to be good when the time-dependence of 
the scale factor is very rapid or, equivalently, 9 is very small as, in this case, the transition 
between sub- and super-horizon regime is very sharp. 

In the following section we compare the approximate value pcst against the numerical 
calculations of /j<^_o, P(t>,2, P<t>A-, Pxfi, Px,'ii PxA- 

4 Numerical results 

For large a = VO/ip and/or small (3 = Hmip the value of pcst significantly underestimates 
the energy density both during and after the inflationary phase, but, in the opposite case, 
the situation improves a lot. In particular, if the relation (2.8) is strictly applied, such that 
60 e-folds requires a = 0.0136, we expect the approximation to be excellent. In figs 3, 4 and 

5 we compare the various definitions of the energy density given in the text for = 3,15 
and 30 respectively. The energy density, defined in units of pmax = (327r^^^)~"'^(3/2)A'^^a^, 
is expressed as a function of the normalized cosmic time T = t/^/89. The grey line is 
the non-commutative energy density which drives inflation, defined in Eq. (3.9). When 
the energy density of the particles generated by the inflationary expansion of the Universe 
becomes of the same order as p^c, it can no longer be neglected, inflation ends and the 
Universe becomes dominated by these massless particles. 

For rapid expansion, the numerical calculation becomes very difficult since the result 
comes from the subtraction of a'^\{vk/a)*\'^ + k'^\vk\'^ and k which, in the case, for example, 
of 60 e-folds, needs more than 100 decimals precision. Hence we have to compute Vk with 
an accuracy better than 100 decimals. As long as the number of e-folds is small, this 
compensation is less severe, but already for = 15 the calculation takes several hours 
(with Mathematica and a precision of 30 digits). 

Let us discuss fig. 3 in some detail. (For this case with Ng = 3, we also show the 
energy spectra in Appendix B.) The energy densities Ptp^n, (purple line, n = 0, blue, n = 2 
and green n = 4) as well always significantly below p^c (grey) for T < 2. 

The estimated energy density (black) underestimates the numerical results for all T. The 
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Piic/pmax 
Pest/Pmax 

P*.o/Pmax 

^ P<^.2/Pmax 

Pi.4/Pnr^>i 

^ P^.o/Pmax 

— Px-l/Pma 

— Px.JPm^:^ 

-2-101234 

T 

Figure 3. The different energy densities for a = 10.6 and A^e — 3. 

— Pnc/Pmax 
Pest/Pmax 

P^.o/Pmax 

P*.2/Pmax 

P4,A/Pm-.K 

^ P^,o/Pmax 

P^.2/Pmax 

^ PxAlPmm 

-2-101234 

T 

Figure 4. The different energy densities for a = 10.6 and A'e = 15. 

gaps in the curves (or sharp raises/decays) correspond to the time intervals when the energy 
density is negative (or to times where it becomes positive/negative). Both, p^^2 ^-nd /30^4 are 
negative when the time evolution is very rapid, i.e. around T = which is very unphysical. 
Note also that at early times, T < —0.5, p0^4 becomes larger than p^^Oi another fact that 
shows how the subtraction scheme cannot be trusted in this regime. The maximum value 
P(^^4 at r < is dominated by the unphysical logarithmic infrared singularity introduced 
by the 4th order subtraction (see spectra of P(f)^4 and in Appendix A). 

For A^e = 15 and Ng = 30 the situation is similar, just a bit more extreme, see figs. 4 
and 5. Now />0 2 remains negative until even later and the infrared divergence of ^,^4 
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Figure 5. The different energy densities for a = 10.6 and A^e = 30. 

becomes more prominent for Ng = 15. For A'^e = 30, it is even much larger than pnc at early 
times, T < — 1 when particle creation should still be very mild since expansion is still slow. 
For both cases p^pfi is well approximated by pest after T ~ — 1. In both cases, all curves 
converge to the estimated energy density at T > 1.5 for the and at T > 3.5 for Px,n- 
Therefore for > 15 the final energy density is well approximated by pest- 

In terms of the dimensionless time T, the time interval of rapid evolution, i.e. inflation, 
always corresponds roughly to —1.2 < T < 1.2. Nevertheless, it is very difficult to calculate 
the energy density for large A'^e since it is hard to numerically obtain Vk with the necessary 
very high precision which is given roughly by (am,/a(l))^ ~ exp(2A^e))) if we consider that 
a^p^O"^ is of order unity (see remark above). Nevertheless, all energy densities converge 
to the same value at late time where time evolution of the scale factor becomes slow. We 
find that this final energy density as a function of the non-commutative energy density is 
approximately given by 



p{t ^ oo) ~ pUT = 2.5) . 



(4.1) 



Note, however, that here we have computed the energy density of the created particles 
from one particle species, while we typically need more than one species to have sufficient 
inflation, depending on the relation of (5 and a. With /3 = \/iV /{a^y/Vhr) we can also 
rewrite eq. (2.10) as 



2N 1 
3 a 



(4.2) 



With a = 10.6, as in our case the true number of degrees of freedom is = (ISA'^g)^. 
The energy density of the created particles, in principle, has to be multiplied by this number, 
which raises p^^ above /?nc at the maximum, T = in the case = 15 and even earlier for 
Nf^ = 30, but e.g. p0_2 remains below p^c until T ~ 1.5 in both cases. Assuming that the 
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true result would lie between p^p^Q and /9<^_2 , as a rule of thumb we can assume that radiation 
domination sets in at about T ~ 1 and therefore back-reaction reduces the effective number 
of e-folds by about 10% (since Erf(l) ~ 0.8). 




-2-101234 



T 

Figure 6. The non-commutative and the estimated energy densities for a — 10.6 and — 60. 



Even if we cannot present precise numerical results for an example with 60 e-folds of 
inflation, we understand what happens qualitatively: for not very much fine tuned values of 
a and /? we can obtain around 60 or more e-folds of inflation if we neglect particle creation. 
Including particle creation, this number will be somewhat reduced but it is difficult to say 
by how much, since all our expressions for the energy density cannot be trusted in the very 
rapidly time-dependent phase. Nevertheless, the created particles, eventually dominate the 
energy density and terminate inflation. This model therefore has a natural, if not very 
graceful, exit from inflation. 

As an illustrative example, we have plotted the background energy density pnc and the 
estimated energy density p^gt for a = 10.6 and = 60 in fig. (6). The black curve has to 
be lifted by a factor A'^ ~ 6 x 10^ in order to account for this number of degrees of freedom. 
This would lift it above the background density already before T ~ and therefore reduce 
the number of e-folds by a factor 2. However, since our calculation of the energy density in 
quantum particle production cannot be trusted in this regime, and since the true maximum 
amplitude of the non-commutative energy density is not very well determined, this result is 
uncertain. If the truly generated energy density remains below /?nc7 

up to T > 1 the model 

can provide close to 60 e-folds of inflation. 

Nevertheless, these estimates tell us that this model has a natural maximum for the 
number of e-folds that can be achieved and which is probably somewhere around ~ 50 
- 60. 



-14- 



5 Conclusion 



In this paper we have considered a model of non-commutative inflation proposed in Ref. [15]. 
We have found that this model leads to a rapid burst of particle creation which we term 
"explosive particle creation". This burst eventually terminates inflation and it is not certain 
whether sufficient inflation can actually be achieved. 

However, if this is possible, particle creation provides an exit mechanism from inflation. 
Inflation terminates at Tend ~ 1 when the energy density of the generated relativistic 
particles takes over. After this, the Universe is radiation-dominated with an energy density 
given by 



We expect interactions to be efficient so that these particles rapidly thermalize and we end 
up with a hot plasma of elementary particles as in standard cosmology. 

We conclude with a comment on the power spectrum of primordial fluctuations. It can 
be shown that it strongly depends on the adiabatic subtraction scheme. If we subtract up 
to the second adiabatic order, the resulting power spectrum is quite flat, p{k) oc log(A;), 
in the UV regime. Therefore, in addition to providing the background radiation energy 
density at the end of non-commutative inflation, this model might also lead to a nearly 
scale invariant spectrum of large scale cosmological density fluctuations. However, to com- 
pute their amplitude and the spectrum in detail, we should push the calculations in the 
far infrared, where adiabatic subtraction is not justified. How to overcome this difficulty 
remains an open problem for a future project. 
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A Adiabatic subtraction 

In this appendix we briefly summarize the adiabatic subtraction scheme, for more details 
see [26]. We assume that the modes can be written in the WKB form 

Plugging this expression into the equation of motion (3.5), one finds that Qk satisfies 






2Qk 




One then solves this equation iteratively to the order n, which concides with the number 
of derivatives of O^. This procedure is valid only in the adiabatic regime, i. e. for k'^ ^ 



l-L^ ,1-1 ~ (rifc/r^fc)^. The term a/ a is to be considered as order two (see [28]). We denote 
the expansion to order n of fifc and of the mode function vt by ^k,n and 

respectively. 

At zero order, we simply have ^kfi = ^ and v"^'^ = v™'^. Removing this contribution is 
equivalent to removing the Minkowski, zero-point vacuum contribution. To find the second 
order, we insert Q.kfi on the left-hand side of equation (A. 2) and expand in d/{ak'^) <C 1 
yielding 

= K' - ^) ■ 

By repeating the operation, we finally obtain the fourth order expression for 0^, 

, / a a(^)o2-2oo2-2a(3)ad-F2a2a\ 
= 2i?a^ «V ]■ <*-^) 

As said before, these expressions are only valid in the adiabatic regime, although some 
authors argue that the subtraction can be employed even when the expansion is strongly 
time-dependent [27]. 

Instead of solving the mode equation for Vk one also could solve the equivalent equation 
for the amplitude = 1/(20^) which follows from (A. 2). One might think that this is 
numerically less involved since one does not have to keep track of the oscillations of the 
modes Vk, see [32]. However, to achieve the very high accuracy needed for our purposes, we 
have found that the numerical problem is as hard as the one for Vk- 
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B Contribution of modes to the energy density 



The energy density spectra, Pjj^^(dp/dlogK) are shown as functions of the dimensionless 
variables T and K. All of them, except q do have regions where they are negative. 
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Figure 7. Contribution to the energy density as function of K and T for N^ — 'i 
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